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Abstract 

We study the evolution of spatial curvature for thawing class of dark energy models. We examine 
the evolution of the equation of state parameter, w^, as a function of the scale factor a, for the 
case in which the scalar field <f> evolve in nearly flat scalar potential. We show that all such 
models provide the corresponding approximate analytical expressions for ^(O^, fJj.) and tt^(a). 
We present observational constraints on these models. 
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I. INTRODUCTION 



About a decade ago, current measurements of redshift and luminosity-distance relations 
of Type la Supernovae (SNe) l| indicate that the expansion of the universe presents an 
accelerated phase 

a a. 

In fact, the astronomical measurements showed that Type la 
SNe at a redshift of z ~ 0.5 were systematically fainted which could be attributed to an 
acceleration of the universe caused by a non-zero vacuum energy density. As this shows a 
result, that the pressure and the energy density of the universe should violate the strong 
energy condition, p^ + 3p^ > 0, where p^ and p^ are energy density and pressure of some 
exotic, unknown and unclustered matter component, dubbed dark energy [4j (see also Refs. 
5|,|6| for recent reviews). A direct consequence of this, is that the pressure must be negative. 

Various models of dark energy have been proposed so far. Perhaps, the most traditional 
candidate to be considered is a non- vanishing cosmological constant [7|,l8|. Other possibilities 



are quintessence 
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k-essence 
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-|l3|. phantom field 14Hl6|. holographic dark energy 



18(, etc. (see Ref. [19[ for model-independent description of the properties of the dark 



energy and Ref. [20j for possible alternatives). 

The first step toward understanding the property of dark energy is to make clear whether 
it is a simple non-vanishing cosmological constant or its genesis comes from other sources 
which dynamically change in time. It is possible to distinguish between these two possibilities 
by taking into account the evolution of the equation of state parameter defined by = 



In what concern to the dynamical dark energy (or quintessence) its physics is described 
by a scalar field, 0, (quintessence scalar field), with canonical momentum 21|. One of the 
main characteristic of the quintessence field is when it rolls the self interacting potential 
curve. It will provide a negative pressure if the potential curve is quite flat. In this way, the 
quintessence scalar field evolves slowly enough to drive the present cosmic acceleration. 

Since the evolution of the quintessence scalar field my be described by the change of the 
equation of state parameter w^, so that we could distinguish two possible situations: the 



case in which du^/dip < and dcu^/dcj) > 0. 
the later the thawing scenarios, respectively 



he former case is referred as the freezing and 



22j(see also Ref. |23| for details). While the 



observational data up to now are not discriminating in the sense that we could not distinguish 
between a freezing or a thawing phases by the variation of the equation of state parameter, 
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it is expect that will be able to do so with the next decade high-precision astronomical 
observations. 

On the other hand, in what concern to the curvature of the universe, today we do not 
know precisely the geometry of the universe, since we do not know the exact amount of 
matter present in the Universe. Various tests of cosmological models, including spacetime 
geometry, galaxy peculiar velocities, structure formation and very early universe descriptions 
(related to the Guths inflationary universe model 24[]) support a flat universe scenario. 
However, by using the seven-year Wilkinson Microwave Anisotropy Probe (WMAP) data 
combined with measurements of Type la SNe and Baryon Acoustic Oscillations (BAO) in 
the galaxy distribution, it was reported that the value for the curvature density parameter, 
Qk = — 0-0057^g;QQgg (68% CL) represents a preferred model, which is slightly closed |25l. |26|. 

In this paper we would like to study some of the consequences that this slightly curvature 
may have on the evolution of the universe, together with the the situation in which the 
thawing cosmological evolution for the quintessence scalar field is invoked. The outline of 
the paper goes as follow, in section II we present the model to be study. Section III, deals 
with the fundamental field equations which allow then and the dynamical system. Finally, 
in section IV we conclude with our finding. 



II. THE MODEL 



The Friedmann equation in which curvature is taken into account becomes given by 

where the Hubble parameter H = a/a, with a dot representing a derivative with respect to 
the cosmological time, a is the scale factor, and the curvature parameter k = 0, +1, and — 1 
represents flat, closed and open spatial section, respectively. Here, we use units for which 
8irG = 1. The total energy density p is given by p = p^ + p m , where p^ and p m are the 
energy density of dark energy and dark matter, respectively. We will assume that these two 
components are conserve separately, satisfying the continuity equations 

pm + 3Hp m = 0, (2) 

and 

pj, + 3#(p + pj,) = p\ + ZHp^l + w^) = 0, (3) 
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where w$ is the equation of state parameter introduced in the introduction. 

We assume that the dark energy is modelled by a minimally-coupled scalar field <ft, where 
the pressure and density of the scalar field are given by 

P*=y-^), (4) 

and 

respectively. Here, V(4>) represents the effective potential associated to the scalar field. 
In term of the scalar field, Eq.flS]) can be written as 

4> + + = 0. (6) 

Equation (jBJ) indicates that the field rolls down the hill in the scalar potential, V(<f>), but its 
motion is damped by a term proportional to H. 



III. EVOLUTION WITH CURVATURE 



We will assume that the scalar field moves in a nearly flat scalar potential, V(4>), quan- 
titatively expressed asj^lj 

?df) <<X ' (?) 

and 

1 d 2 V 

ia « L (8) 



V 



71) and (El) indicate 



The constraint given by Eq. ([TJ) ensures that ~ —1, meanwhile Eqs.( 
that (l/V)(dV/d(/)) is nearly constant (27)]. In the nomenclature of Ref.[22|, these are thawing 
models, i.e. du^/dcj) > 0. 

From the Friedmann equation we get 



where the density parameters are Q^, = p^/(3H 2 ), Q m = p m /{3H 2 ) and Qk = —k/(aH)' 
respectively. 
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Following, a similar technique developed in Ref. [271]. Eqs.((T]) and ()6]) can be expressed in 
terms of new variables x, y, A, and fl k , defined by 

x = <///VE, (9) 

(10) 

(11) 
(12) 



y = vV(0)/3i/ 2 , 

n k = k/h 2 , 



where K = —ka 2 and a prime denote the derivative with respect to In a, and V<j, = 
dV{<j))/d<j). 

The density parameter is expressed in terms of the variables x 2 and y 2 in such a way 
that 

n^ = x 2 + y 2 , (13) 
while, the equation of state parameter is given by 

7 = 1 + w<f, 



2x 2 



(14) 

x l + y z 

Eqs.((T]) and ([6]) can be written in terms of the new variables Eqs.()9|)- (|T3|) . so that we get 



x 



/3 3 
-y 2 H — x 
T 2 



1,2 2 

1 + x — y 

y 3 



. 3 3 

■ x yr y+ 2 y 



1,2 2 

1 + x — y 

y 3 



where 



y = 

\' = -V6X 2 (T-l)x, 

n' k = n k (i-n k + 3[x 2 - y 2 }) 



\d(f> J 



(15) 

(16) 

(17) 
(18) 

(19) 



In the thawing model we have ~ — 1 and thus the 7 parameter satisfies 7 = 1 + <C 1. 
Therefore, it is useful to express Eqs. (TT5"j) - (1T81) in terms of 7, Q^, Q k and A, respectively. 
We obtain 



i = -3 7 (2 - 7) + A(2 - 7) ^3W^ (20) 

% = 3(1 - 7)^(1 - fy) - n^fifc, (21) 

Q' fe = fi fc (1 - fi fe + 3fi^(7 - 1)) (22) 

A' = -V3A 2 (r-l)V?V (23) 
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At this point we would like to stress two assumptions that we are considering: the first 
is 7 1, which corresponds to ~ — 1, as discussed previously. The second assumption 
we make is that the scalar field begins with an initial value in a potential which is nearly 
flat. In this way, following j^, we assume that A is approximately constant, so that 



An 



-(l/V)(dV/d(f>) 



« 1, 



(24) 



<j)=4>o 



where Ao is a small constant evaluated at = <po, the initial value of the scalar field which 
corresponds to when it stars to roll down the potential. 

Let us first to consider the evolution of the system using initial values for the curvature 
Qk — 0.005. The result is shown in figured] The same graph, but now with a value for the 
curvature Qk — —0.005 we get a similar plot with a small difference when compared with the 
previous case. The area in between the curves expands a continuous range of values of the 
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FIG. 1: Numerical results for 7 as a function of the fractional density parameter f2^, for nearly 
flat potentials. We have used the present values tlko = ±0.005 and fi^o = 0.7, with Ao = 0.01.. 



curvature parameter, well inside the current observational constraints. In order to see this 
situation more clearly, we plot in figure [2] the projection of curves in the f^-fi^ plane. We 
observe that a large region exist, even for small values Qk- As a complement, in figure [3] we 
show the degeneracies in the variation of 7 with respect to the curvature. All these figures 
are the result of a numerical integration of the system of Eqs. ( 12011231) with Ao = 0.01. 

Motivated by the present value measured for the curvature parameter we make the as- 
sumption that the curvature is a small parameter, i.e. <C 1 along of the all story of the 
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FIG. 2: Numerical results for £1% as a function of the fractional density parameter Q^, for nearly 
flat potentials. We have used the present values O^o = ±0.005 and Qm = 0.7, with Aq = 0.01. 



0.005 



cf 0.000 - 



-0.005 



0.05 



0.10 

y 



0.15 



0.20 



FIG. 3: Numerical results for as a function of 7, for nearly flat potentials. We have used the 
present values Q^o = ±0.005 and tl^o = 0.7, with Aq = 0.01. 



evolution of the universe. From Eqs. ( 120]) and (12 ip we can write 

c?7 _ -67 + 2A O a/37^0 AoOfc^/37^ 



(25) 



where we have expanded and maintained the lowest order terms in and 7. Taking the 



boundary value 7 = at fi^, = (see Ref.[27|). The resulting solution is 
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1296 



12 + 



fi<?i 2(Q^ — 1) y fi 9 



F(D, k , fi^) 



where 



-F(fi fc , fl q 



\1) (12 + ^)tanh- 1 v ^ 



(26) 



(27) 



Note that in the limit of a flat universe, i.e., fife — > 0, we recover the expression given in Ref. 



3- 



In the same way we can derive an approximate solution from Eqs. fT2T]) and (1221) under 
the same approximations (7 -C 1 and fife <C 1) 

(2fi fe fi fe (l - 3fi,i) 



rffij, 30 1 — fi„ 



from which we get 



fife — fifco 
From Eqs. ([26]) and ([29]) we get 



" 1 - v 


2/3 




. 1 — fi</>o . 


1 





1/3 



(28) 



(29) 



\2 

1296 



12 ^fco 





2/3 








/ 



(fi + 1; 



2(^-1)^ 



G(fi</,) 



(30) 



where the function (7(0^) is given by 



fit., 



A'O 



1-n 



<0() 



2/3 ' : 



tanh 1 \/fL 



/ 



We can use equation ( 121]) to solve for fi^, as a function of a and thus determine w^a). 
Taking the limit 7 1 and fife *C 1 in equation f ]2~T]) gives the following solution 



= [1 + (fi^ 1 - 1) a~ 3 ] 



31-1 



(31) 



where fi^o and fifeo are the present values of fi^,, fife, respectively, and we take a = 1 at 
present time. Combining Eq. f[3"Tj) with Eq. ( 12U|) we obtain an approximated solution for 
fife (a). Then, with the explicit expressions for fi^(a) and fife (a), and by using Eq. ( 126]) we 
get explicitly the equation of state parameter u^, as a function of the scale factor a, i.e. 
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In Fig. @ we show the dependence of the parameter as a function of the scale factor 
a, for different values of the curvature parameter Qk with wn = —0.95 and f2^o = 0.7. Note 
that w^{a) is not sensible to the value of Qk — (see Ref. 27|). 

We should mention that if we look for numerical solution to our set of dynamical Eqs., 
in which a scalar potential, such that V{4>) ~ <f) 2 , <f)~ 2 , exp[— 0], etc, is used, we observe 



that there is no much changes when them are compared with that shown in Ref. 27J, where 
Slfc = was taken into account. 
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FIG. 4: Our anality results for the evolution of the parameter as a function of the scale factor 
a for nearly flat potentials. Here, we have taken the values wq = —0.95, Q^o = 0.7, Qko = +0.005 
(dot line), fi^o = (solid line) and £1^0 = —0.005 (dash line), respectively. 



Having an approximated expression for iu^(a) we can use it to perform a Bayesian analysis 
using SNIa observations, BAO distances and CMB shift parameter. In this work, we use the 
Supernova Cosmology Project Union sample 28j, having 307 SN distributed over the range 



0.015 < z < 1.551. We fit the (theoretical) distance modulus fJ,(z) t h defined by 

H d L (z) 



n(z) th = 5 log 



10 



+ Mo, 



(32) 



to the observational ones n(z) b s - Here Hq = 100/ikm s 1 Mpc 1 is the Hubble constant and 
the luminosity distance is defined by di(z) = (1 + z)r(z) where 

dz' 



r\z) 



:Sinn-\/|Qs 



(33) 



and /io = 42.38 — 51og 10 /i. Sinn(x) = sin re, x, sinh x for < 0, Qk = 0, and Qk > 
respectively. The second major input for parameter determination comes from the baryon 
acoustic oscillations (BAO) detected by Eisenstein et al. |29[. In our work, we add the 



following term to the x 2 of the model: 



Xbao 



(A- A 



BAO, 



OA 



(34) 



where A is a distance parameter defined by 



.4 



czbao 



2, \ CZBAO 

r (zbao) 



1 V3 



(35) 



H{z BAO ) 

and Abao = 0.469, a a = 0.017, and zbao = 0.35. The CMB shift parameter R is given by 
30] 

R{z m ) = 



(36) 



Here the redshift (the decoupling epoch of photons) is obtained by using the fitting 
function [31 ] 

z, = 1048 [1 + 0.00124(fi 6 /i 2 )- - 738 ] [1 + gi (n m h 2 ) 92 ] , (37) 
where the functions g\ and g 2 are given as 



9l = 0.0783(^ 2 )-°- 238 (l + 39.5(^ 2 ) - 763 )" 1 
g 2 = 0.560 (l + 21. l(fi 6 /i 2 ) x - 81 ) -1 . 



(38) 
(39) 



The WMAP-7 year CMB data alone yields R{z 
sponding Xcmb as 



XCMB 



-- 1.726 ±0.018 
(R(z*) - 1.726) 2 



321 ] . Defining the corre- 



(40) 



0.018 2 

one can deduce constraints on f2^ , an d ^fco- A joint analysis using SN+BAO+CMB 
leads to the best fit values showed in FigJH where we see the cross section of the x 2 function 
in terms of the parameters Q^o, ojq and Qko- The two horizontal lines indicate the 90% and 
99% confidence range for each parameter. 

The analysis shows that considering thawing quintessence with an explicit curvature 
term is consistent with observations. This is exactly the conclusion of 27] for the flat case 



in quintessence. However, as was demonstrated in [33j, relaxing the slow-roll assumption, 
the equation of state parameter for different thawing potentials looks appreciably different. 
In the following, we consider both quintessence and Tachyon field models, and two scalar 
field potentials; V = <fi and V = <p~ 2 . In figure ?? we show the integration of the field 



equations for current values of the curvature Vt 



k0 



±0.006 and ft 



<;>0 



0.72 for all the 



models. The potential are characterized by T = and T = 3/2 respectively, along the initial 
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FIG. 5: This panel shows the \ 2 function computed using the approximate analytical solution of 
w^(a). The best fit parameters for O^Qi w o an d Qko- The upper dashed line indicate the 99% 
confidence range for each parameter and the continuous line below indicate the 90% of confidence 
range for each parameter. 

conditions A™, ~ 1 (assuming that the potential is not flat) and 7i„j ~ (the equation of 
state parameter can vary from its freezing state (w — — 1) until today.) 

IV. CONCLUSIONS 

In the present work we have studied the thawing dark energy scenarios in which the effect 
of curvature was taking into account. 
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FIG. 6: Here we show the numerical integration for two scalar field models; quintessence and 



explicitly Q^o = ±0.006 with a current value Sl^ = 0.72. The upper two curves correspond to a 
Tachyonic model with potential V = (p. Although the values of the curvature parameter are very 
small, the separation of the curves increases with evolution and are appreciable here. The next 
couple of curves correspond to a quintessence scalar field model with potential V = <fi. The third 
set of two curves (which are closer each other than the previous ones) correspond to a Tachyonic 
model with potential V = 4>~ 2 - The bottom two curves correspond to a quintessence model with 
V = (j)~ 2 . All these models have = 0.72. 

We have plotted numerically trajectories in the (7, fi^), (flk, ^0) and 7) for a 
potential nearly flat. 

We have shown that all such models converge to a common behavior and we have find 
the corresponding approximate analytical expressions for 7(^0) given by Eq.( l30|) and for 
Wffyia) in the cases when 7 1 and fij. < 1. Here, we noted that an analitical solution for 
w^ia) is not very perceptible to the value of Vt k 7^ 0. A Bayesian analysis using SNIa data 
was performed to constraint the best fit parameters using our analytic function, w^{a). This 
analysis shows that current data does not rule out the model. In this way, the motivation is 
to see whether one can distinguish thawing dark energy models from Vtk 7^ models using 
this method. 



Tachyonic. We consider two scalar field potentials V = <f> and V = <fi 



2 and in each case we use 
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